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Abstract 

We study the quantization of many-body systems in three dimensions in rotating coordinate frames 
using a gauge invariant formulation of the dynamics. We consider reference frames defined by linear 
gauge conditions, and discuss their Gribov ambiguities and commutator algebra. We construct the 
momentum operators, inner-product and Hamiltonian in those gauges, for systems with and without 
translation invariance. The analogy with the quantization of non-Abelian Yang-Mills theories in non- 
covariant gauges is emphasized. Our results are applied to quasi-rigid systems in the Eckart frame. 

1 Introduction 

The problem of quantizing a many-body mechanical system in a rotating reference frame is of interest both 
by itself and for its possible applications to specific problems in, e.g., molecular and nuclear physics. In this 
paper we study the quantization of many-body sytems in three dimensions in rotating coordinate frames, 
using a gauge- invariant formulation. The two-dimensional case was considered in a previous paper pQ, in 
which the method was developed in detail and a close parallel with the quantization of electrodynamics in 
non-covariant gauges established. The main lines of the method are the same in both cases. Due to the 
non-Abelianity of the rotation group in three dimensions, however, the technical treatment of the systems 
considered here is considerably different from the planar case, the differences being already apparent at the 
Lagrangian level as discussed in the following section. 

We consider systems of N spinless particles interacting through central potentials. Since the underlying 
dynamics are rotationally symmetric, the coordinate transformation from a space-fixed reference frame to a 
rotating one with the same origin is a time-dependent symmetry transformation, or gauge transformation. If 
the dynamics are described in terms of a gauge- invariant action, since we know how to quantize a mechanical 
system in a space-fixed coordinate frame, we can perform a gauge transformation in order to obtain the 
quantum theory in a rotating frame. Gauge invariance implies that both theories are physically equivalent. 

Rotating frames are often defined implicitly, by restrictions on the trajectories of the system in that frame. 
In the gauge- invariant approach to the quantization in rotating frames, such restrictions are incorporated 
into the theory as gauge conditions. The action is then given in terms of degrees of freedom that are not 
independent, but must satisfy certain functional relations. This situation is familiar from the theory of 
gauge fields 0EI], where the vector potential A(t, x) may be required to satisfy such relations as V • A = 
(Coulomb gauge), or n- A = (axial gauge), at all times t. In this paper we consider only gauge conditions 
depending linearly on the particles coordinates, which are most useful in practical applications involving 
perturbative expansions. We do not consider quadratic gauge conditions, in particular, because we expect 
the formalism in those gauges to be considerably more complicated, in view of the results of [P in the simpler 
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two-dimensional, Abelian case. Furthermore, the quadratic gauge conditions most common in the literature 
0J|S] are those defining the instantaneous principal axes frame, in which the total angular momentum of the 
system is strongly coupled to the other degrees of freedom through the inertia tensor and the Coriolis terms. 
Without a strong physical motivation for quadratic gauges, we have no reason to pursue that technically 
more involved approach here. As discussed in PP, however, there is no problem of principle to deal with 
those and other kinds of gauge conditions within the formalism espoused in this paper. 

We closely follow the approach to non- Abelian Yang-Mills theories in non-covariant gauges of EJ, 
stressing throughout the paper the strict formal similarity between our results for many-body systems and 
the corresponding ones in jS] for Yang-Mills theories. Our goals are both to illustrate the formalism of gauge 
theories in the more familiar context of mechanical systems, and to apply the gauge-theoretical techniques to 
the quantization of three-dimensional TV-body sytems in rotating frames. Previous treatments of the latter 
problem within a gauge-invariant approach have been given in 0] and references therein. A gauge theory 
of rotations and internal motions of deformable bodies, including classical and quantum A-body systems, is 
developed in \7\ (see also [5]) from a point of view different from ours. Non-gauge-invariant treatments can 
be found, e.g., in |S| in the context of nuclear physics, and in |51 110| in molecular physics. 

The outline of the paper is as follows. In section [2] we describe the class of systems to be considered 
throughout the paper, and their formulation in terms of a Lagrangian invariant under time-dependent 
rotations. Their quantization in a space-fixed frame is given, and shown to be equivalent to the non- 
gauge-invariant formulation. The central results of the paper are given in section |3 where we consider the 
quantization in rotating frames defined by linear gauge conditions. We discuss in detail the commutator 
algebra for both linear and angular momentum operators, and give explicit realizations of that algebra in 
terms of differential operators. Those operators are used to construct the Hamiltonian in terms of position 
vectors referred to the rotating frame, and their conjugate momenta. The elimination of orientational 
degrees of freedom from the formalism is subsequently carried out, and the resulting Hamiltonian and its 
Weyl-ordered form and related quantum potential are obtained. As emphasized throughout, by describing 
a many-body system from a rotating frame defined by gauge conditions, we are introducing orthogonal 
curvilinear coordinates in configuration space. The singularities of those coordinates occur at the Gribov 
horizons where the gauge conditions become ambiguous. Gribov ambiguities [111 l?fl Efl IT] are discussed in 
detail in relation to the construction of the inner product in the reduced state-space of the system. 

In sectionQ]we extend the results of section[3]to translation-invariant systems. We show how the gauge- 
invariant approach can be used to describe a mechanical system in a reference frame in an arbitrary state of 
rotation and translation. In particular, we obtain explicit results for the quantization of A-body systems in 
rotating frames with origin at the center of mass. Those results are then applied to quasi- rigid systems in the 
Eckart frame ^1 EID in section where two simple three- and four-body examples are briefly discussed. 
In section we give our final remarks. Some complementary material is gathered in the appendices. 



2 iV-particle system 

We consider a system of N spinless particles with central interactions in three dimensions, described by the 
Lagrangian 

1 N 

£ = C N + C T t , £n = 2 X! m ^a - V 

° =1 (1) 
N N X 2 X 2 

V= V ocp(\r a -r p \) + J2u(r a ) , £ rt = -(eAeJ = - e , e-e = l. 

a<0=l a=l 

The potential energy V is chosen for concreteness to include only one- and two-body interactions. Our results 
do not depend on that fact and apply equally well to more general rotationally invariant potentials. If the 
one-body potential U = 0, C is invariant under the group of Euclidean motions of three-dimensional space. 
In this and the following sections we consider U ^ and focus on the non- Abelian group of three-dimensional 
rotations, deferring the discussion of translation invariance until section^] Besides the kinetic and potential 
energy for the A-particle system C also contains the Lagrangian £ rt for a free rigid rotator, described by 
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a unit vector e. This rotator is not coupled to the particle system, so it does not affect its dynamical 
evolution. £ rt can be made to vanish by letting the rotator's inertia moment X — > oo while keeping constant 
the magnitude of its angular momentum s = Xe A e. Whereas C T t does not play any role in C as given in 
(|T|) , it will serve as a source of angular momentum for the particle system in the gauge-invariant formulation 
to which we now turn. 

C is invariant under time-independent rotations of the coordinate frame. In order to make C invariant 
under changes of arbitrarily rotating coordinate frames we apply the Yang-Mills construction JSj to (QJ. 
We introduce a 3 x 3 real antisymmetric matrix thus adding three new degrees of freedom to the system, 
and postulate the following transformation laws under rotations of the coordinate frame, 

r' a = Ur a , e ' = Ue , £' = U£U* + UU** , (2) 

with U a time-dependent real orthogonal matrix and U 1 its transpose. These are the gauge transformations 
of the system. We define the covariant time-derivative D t r a = r a — £r a , which transforms like a vector 
under gauge transformations, (D t r a )' — U(D t r a ). Analogously, D t e = e — £e. We can, equivalently, use 
instead of the matrix £ the axial vector £ dual to ^ = £ijk£,j, whose gauge transformations can be 
derived from (J2J). I n terms of £ covariant derivatives take the form D t r a = r a — £ A r a . Substituting time 
derivatives in Q by covariant derivatives we obtain a Lagrangian invariant under time-dependent rotations 
of the coordinate frame. Explicitly, we write, 



C — Cn + £rt j 

1 N i N i N ~ 2 „ N 

£n = ^ X! m a {D t r a ) 2 - V = - ^ m a r 2 a + - m a A ~ £ ■ m a (r a A r a ) - V 



= 1 a=l 



(3) 



Crt = \ (e A (D t e)f = - (D t ef = ~£ + - (t A e)* - T~t ■ ( e c 



where the potential energy V is defined in (Q. £ is exactly invariant under the gauge transformations (J2J 
and, in fact, Cn and C r t are separately invariant under Q. 1 

Cn in has the form of a Lagrangian for an ./V-particle system described from a coordinate frame 
rotating with angular velocity — £ with respect to the laboratory frame |14) . Notice, however, that £ is 
a dynamical variable describing the coupling of the particles and the rotator to the inertial forces. The 
equations of motion for r a , e, and £ derived from C are 



m a D t D t r a + V Q V = m a r a — 2m a £ A r a — m a £ Ar„- m Q £ A yr a A £j + V Q V = , (4a) 

D t D t e + (D t ef e = with ee = l, (4b) 
dC N 

= = V m a r a A (D t r a ) + Xe A (D t e) = . (4c) 

d € a =i 

In the equation of motion l|4a|l for r a the terms due to the Coriolis, azimuthal and centrifugal forces ^1 are 
apparent. As a consequence of the rotational invariance of C the total angular momentum of the system is 
conserved, dj/dt — with, 

N 

j = l + s , l = ^2 maVa A iPt r a) i s = Xe A (D t e) , (5) 

a=l 

Clearly, the vector j can be time-independent in every rotating reference frame only if it vanishes. This 
is expressed by (|4cf) . which can be rewritten as j — 0. Since in general C is not invariant under separate 
rotations of {r a } and e, I and s are not separately conserved. Rather, from l|4a|l and i|4b|l they are seen to 
be covariantly conserved, 

D t l = , D t s = , (6) 

1 Notice that, unlike the two-dimensional (Abelian) case [3, we cannot add external-source terms to C without breaking 
gauge invariance, so we have to incorporate the source into the theory as a dynamical degree of freedom. That is the motivation 
for including £ r t in C We stress here that £ r t is gauge invariant and therefore it is not a gauge fixing term. 
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with D t l = 1 — £ A I. From ©, the magnitudes of I and s are conserved and frame- independent, but their 
directions in space are time-dependent. Only in the lab frame (in which £ = and D t = d/dt, as discussed 
below) are I and s conserved. 

Since the system is gauge invariant we can fix the gauge by imposing a set of conditions of the form 2 
© a ({r Q }, £, e) = 0, a = 1, 2, 3, which is equivalent to selecting a rotating frame in which the trajectory of the 
system ({r a (t)}, £(i), e(t)) in configuration space is constrained by the relations <3 a ({r a (t)}, £(t), e(t)) = 0. 
The functions 3 <3 a , a = 1,2,3, can be chosen arbitrarily, as long as any trajectory ({r' a },£\ e') can be 
transformed into a new one ({r Q },^, e) satisfying © a = 0. The new trajectory must be unique, in the sense 
that no other trajectory obtained from ({r^}, e ') by a gauge transformation satisfies the gauge conditions. 
Otherwise, the gauge is said to be ambiguous Supplementary conditions must then be imposed to fix 
the ambiguity. 

2.1 The laboratory frame 

Given any trajectory of the system ({r a (t)}, e(t)) by means of a gauge transformation we can obtain 
a physically equivalent trajectory with = U^U^ + UU^ = 0. Indeed, given the antisymmetric matrix- 
valued function of time £{t), there is always an orthogonal matrix U(t) satisfying U^U = — £. The condition 
£ = is then admissible as a choice of gauge for the system, which corresponds to selecting a non-rotating 
coordinate frame referred to as the "laboratory frame." 

We denote dynamical quantities in the laboratory frame by lower-case symbols, except for the Lagrangian 
and Hamiltonian. In this gauge the Lagrangian J2J reduces to £ is invariant under separate rotations 
of {r a } and e, leading to the separate conservation of the angular momenta I = Yla=i m a r a A r a of the 
system of particles and s — Te A e of the rigid rotator. The equation of motion ljlc|l for £, which cannot be 
obtained from 1|T|1. must be imposed on the system as a constraint 2 , j = I + s = 0. In the Hamiltonian 
formulation in this gauge, this is a primary first-class constraint |16| . not leading to further secondary ones. 

The quantization in the gauge £ = is canonical. In units such that h — 1 we have, 

N 1 1 

h = h n + Hn, n N = Y, + V . ft rt = ^ . 

a— 1 a 

[rauPpj] = i5 a p&i 3 , P a = [sj, Sj] = ie i]k s k , (7) 

N 

{<M)= / IpV^* ({r Q },e)V({r Q },e) , 

J (3=1 

with the commutators among r a and p a not shown in J7J) all vanishing. The first-class constraint is imposed 
on the state space \n^, j\tp) = 0. Since both I and s are constants of motion, this constraint is clearly 
consistent with the dynamics. We see that the quantized theory in the £ = gauge is completely analogous 
to Yang-Mills theories in the temporal gauge [HI El El- The constraint fixing the value of j, in particular, 
is the equivalent of the non-Abelian Gauss law. In the constraint equation s plays the same role as the 
fcrmion color current in Gauss law. 

If in J2J we let X — > oo with s 2 fixed, 7i r t — ► and the rigid rotator drops from the Hamiltonian, entering 
the dynamics only through the constant value of s in the constraint. Thus, JZJ) describes in that limit an 
iV-body system with interaction potential V in the sector of fixed angular momentum I = — s (i.e., the null 
eigenspace of (I + s) 2 ). Due to gauge invariance, the same must be true in any other gauge. 

2 The letters a, b, c, d are used to index quantities which are not necessarily tensor components, such as <3 a - Summation over 
those indices and their ranges of variation are always explicitly indicated. We only use the summation convention for tensor 
indices, which are denoted by latin letters i,j,k,l, ■ ■ ■ and always run from 1 to 3. 

3 In general, (S a are functionals of the trajectory ({r a (t)}, £(t), e(t)). 
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3 Linear gauge conditions 



In order to fix a reference frame we need to impose three gauge conditions. The simplest gauge conditions 
involving the coordinates of the particles depend linearly on {r a }, and do not involve £ or e. As discussed 
below, linear gauges are relevant in the context of perturbative or semiclassical expansions. The general 
form of the linear gauge conditions is, 

N 

&a({r a }) = ^2 m a T aaj r aj = , a =1,2, 3, (8) 

a=l 

with T aa j a set of 9 A constants defining the functions & a . We denote dynamical quantities in this gauge by 
capital letters, in particular the position vectors R a , E and the angular momenta L and S, as opposed to 
the corresponding quantities in the gauge £ = (the laboratory frame) which are denoted r a , e, I, s. Thus 
& a ({R a }) = but, in general, 6 a ({r Q }) ^ 0. The gauge conditions © select a reference frame rotating so 
that the linear combinations of coordinates & a vanish for all t. If we choose, for instance, all coefficients in 
||SJ| vanishing except for Thy — ^i\z = ^32Y = 1, the coordinate frame must rotate together with particles 
1 and 2 so that 1 is on the X axis and 2 on the X — Z plane for all t. The formalism in these linear gauges 
is entirely analogous to that of non-Abelian Yang-Mills theories in linear non-covariant gauges, such as the 
Coulomb or axial gauges, in which the fields are also constrained by linear relations [Hj(see also pfl 1171 IT%] V 
For the functions & a ({r a }) to be admissible as gauge conditions they must not be rotationally invariant. 
The variation of & a under an infinitesimal rotation is S& a = £l a kS0k, with 

JV 

QaidRy}) = ^ m P^af3j£jikRf3k , a = 1, 2, 3. (9) 

/3=1 

The requirement that & a must not be invariant under infinitesimal rotations is therefore satisfied if the 
matrix Q a i is not singular on the gauge manifold. Thus, the following equations must be simultaneously 
satisfied, 

& a ({R a }) = 0, a = l,2,3, det(Q 6i {{Ra)}) ^0 , (10) 

except possibly at exceptional configurations at which the gauge is singular, det £3 = 0, such as R a = for all 
a. Furthermore, without loss of generality, we assume that the gauge coefficients have been orthogonalized 
so that 

N N 

^ m a r aaj T baj = 9\ 2 a S ab , ^1 = 22 m a T aaj T aa j > , 1 < a, b < 3 . (11) 

ol—1 a — 1 

The gauge transformation from the gauge £ = to the gauge & a = is of the form 1(3}, 

R a = Ur a , E = Ue , £ = UU 1 . (12) 

The orthogonal matrix U is parametrized by three angles {# a }a=i- Although our approach and results do not 
depend on any specific parametrization of the rotation group, some parametrization-dependent quantities, 
such as the momenta pe a conjugate to 9 a , are physically meaningful and play an important role in some 
intermediate calculations. All the information we will need about the parametrization of U is encoded in 
the matrices A and A defined by 

|^C/t =Aa . T . 5 U ^ =Xai T i , a = 1,2,3, (13) 

where the Tj are the standard generators of the so(3) algebra, (Tj)ik — £ijk- The three matrices dU /d6 a U\ 
a = 1, 2, 3, must be a basis of so(3) for all values of {#b}, if the parametrization is to be well defined. Thus, 
the matrix A a j is invertible and, analogously, so is X a i- From the unimodularity of U it follows that 
U^TiU = UijTj and then, from (|13fl . \ a j — A a iUij. We can express £ in terms of 6 a and their time 
derivatives as 

3 3 

&fc = ^ (taKjSjjk , or gj = ^ 9 a A aj , (14) 

a— 1 a— 1 
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£ in this gauge can also be written in terms of R a , E, and their time derivatives from the constraint equation 
L + S = 0. The resulting expression, unlike l|14|) . would be valid only within the constrained subspace. 

Through (|12fl . the gauge conditions determine the time dependence of {9 a } so that, given a trajectory 
({r a (t)},e(t}) of the system in the gauge £ = 0, we have & a ({R a (t)}) = 6 a ({U({9 a (t)})r a (t)}) = for 
all t. We view (|12ll as a coordinate transformation in configuration space, specifying the new coordinates 
{R a ({r a })}, {9 a ({r a })}, E({r a },e) in terms of the original ones {r a },e. The number of independent 
variables is the same in both sets, since the N position vectors {R a } are restricted by the three linear 
conditions & a {{R a }) = 0. From 1(12(1 we then have, 

9R/3i r jj . dU ik , s 

Substituting l|15|l into the relation d& a ({R/3})/dr a j = 0, and using the definition 10 for 0, the assumption 
(f 1 01) that it is invertible on the gauge manifold, and the antisymmetry of (dUik / dr a jUik) in i and I, we 
obtain the relation 

-5— —Ulk = ^2 £ UmQm 1 a m a^aanU n j , (16) 
UT a 7 -, 
J a—1 

which expresses dU/dr aj U^ in terms of {Ry} and {8 b }. This expression characterizes the dependence of 
U on {r a }, and will be important below, especially in the discussion of angular momentum (see sect. I3.l|) . 
Unlike the two-dimensional case pQ in which the explicit form of U is easily found, 1)16(1 does not give us 
explicit information about possible Gribov ambiguities of this gauge. Those ambiguities are analyzed below 
(sect. l!0)l . in connection with the derivation of the Hilbert-space inner product in this gauge. 

The Lagrangian in this gauge is given by C in Q with r a and e substituted by R a and E, according 
to our convention. Due to the relation (|14f) between £ and 9 a , we can use {R a }, {9a}, E as dynamical 
variables, the Lagrangian in terms of them being obtained by substituting £ = 17Z7' in Formulating the 
theory in those variables, however, would result in momenta pe a conjugate to 9 a which are linearly related 
to J = L + S, not just L. Furthermore, in the Hamiltonian formulation we have [Li, E] ^ = [Ji,E]. 
Thus, the gauge transformation (|12f) "mixes" the particle degrees of freedom {R a } and {9 a } with the rotator 
degrees of freedom E. We can avoid such mixing by describing the rigid rotator in terms of its position 
versor in the lab frame. Once the dynamical variables in have been appropriately capitalized, we set 
£ = UU' and E = Ue to obtain, 

X ■ 2 

C = C N +£ rt , C It = -e , (17) 

N 2 N 3 N 3 

Cn = — 2_j m aR a + 2 m " (^a^ij ~ RaiRaj) ^ A ci Adj9 c 9d ~ m a eijkR a j Rak A ci 9 c — V . 

a — 1 a — 1 c,d—l a—1 c— 1 

This Lagrangian, with the gauge conditions (JHJl holding as strong (operator) equalities and the constraint 
J = valid as a weak (state space) equality describes the same dynamics as ©. The formulation based on 
(|17fl . with {R a }, {9 a }, e as dynamical variables, closely follows the treatment of non-Abelian Yang-Mills 
theories in non-covariant linear gauges given in [S]. 

3.1 Angular and linear momenta 

In the quantum theory in the gauge £ = 0, as discussed in section \l. II the angular momentum operator I 
satisfies the usual commutator algebra. Using 116|) . the definition © of £2, and the unimodularity of U, we 
obtain, 

N qjj 

[h, Ujk] = 2^ £ ilm T al - q J = i£iknUj n . (18) 
a—1 

From (JTSJ, using L = Ul and R a = Ur a , we get, 

[Li,Ujk\ — i^ijnUnk , RoLj\ l-^ii Rotj\ , (19) 

[Zt,Zj] = iSijklk , [h, Lj] = , [Li,Lj] = —ieijkLk . 



6 



As expected, the particle position vectors R a in this gauge are rotation invariant. The commutators among 
components of I and L are the same as for a rigid body, with I the angular momentum in the laboratory 
and L in the body frame. Furthermore, taking into account the commutators Q for s, [s, U] = 0, S = Us, 
and J = L + S = Uj , we have, 

[Si, Sj] = iEijkSk , [Li, Sj\ = —iEijkSk , (20) 

Notice that [L, S] ^ 0, due to the dependence of S on the angles {0 a }. The classical expressions for L and 
pe a follow immediately from the Lagrangian (|17J) and 

N N 

P9a = —^aiLi , Li = m a EijkRajR a k ~ J]] ™q (R a ^ij ~ RaiRaj) £j , (21) 

a— 1 a—1 

with | given by <H3Jt. By using the identity J2l=i dA ci /d9 a c = A ai + K aj E ijk Y^x^cKk, which follows 
from the definition of A ai , the equation of motion for 9 a from the Lagrangian l|17|l can be reduced to 
the form D t L = 0, in agreement with ©. In the quantum theory, 

a— 1 

Similarly, pg a — —X a ih and U = J2 a =i Ka^/^a- Eqs. dt - jHI are completely analogous to the relations 
among color currents in Yang-Mills theories in non-covariant gauges (see eqs. (4.44)-(4.48) and (4.55) in 0). 

In the classical theory we obtain the momenta P a conjugate to R a by differentiating l|17|l (or, equiva- 
lently, J3J)) with respect to R a under the constraints & a ({R a }) — & a {{R a }) = to obtain, 



_ / 3 i \ 

Poii WlaRai THafij f ^-ijkRoik ^ ^ ^2 ^bj^bai 1 (^^) 

with Qbj and 9^ defined in © and These momenta are consistent with the gauge condition, since 

they satisfy 

N 

= Yl r «« P « = & - ({ p a /m a }) . (24) 
0=1 

From the transformation l|12l) we can obtain the relation between the velocities {r Q } in the gauge £ = 0, 
and those in the gauge & a — 0, {Ra}, {da}- Correspondingly, we can express the momenta {p a } in one 
gauge in terms of the momenta {P a } and L in the other, 

/ 3 \ N 

Paj = U kj P ak + ^ m a r aak £l^(L n - A n ) , with A„ = ^ ^n Vq R lv P iq ■ (25) 

\ a=l J 7=1 

The quantity 4 A„ defined by this equation has the appearance of an angular momentum but, as shown below, 
it does not satisfy the so(3) commutation relations in general. With the transformation (|25|l for momenta 
we obtain from 7ijv in Q) the classical Hamiltonian for the particle system in this gauge, 

N 3 

n N = Y. ^T p l + \ E Klttr a ^{ Li A^ A,) + V . (26) 

a—1 a—1 

The Hamiltonian Tt It for the rigid rotator is clearly the same as in 0. 



4 There should be no possibility of confusion between the 3-component operator A defined in 1251 and the 3x3 matrix A a i 
defined in 1131 . 
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The transformations l|12fl and (|25[1 can be inverted, to express P a , L, R a , and 9 a in terms of p a and r a . 
Using the Poisson brackets J7J, we then get the Poisson brackets in this gauge. Alternatively, they can be 
found as Dirac brackets JHI relative to the set of second-class constraints 6 a ({R a }) = = & a ({P a /m a }), 
a = 1,2,3. The results, written in the notation of quantum commutators, are 

[R a i, P/3j] = i I S a pSij — -zJ^F aai T a pj J . (27) 

V a=l / 

All other commutators among R a and Pp vanish, and [L, R a ] = = [L, Pp]. From l|27|l we get, 

[6 a ({i? Q }), P Pj ] = = [e a {{P a /m a }), R 0j ] = [An, & a {{Ra})} , (28a) 

N 3 . 

[Aj, Aj] — iCijkA-k 2 ^ ^ ^ ^ ^y>2 ^fla»i(^™iiQaj ^jrnnQai )Pan ■ (28b) 
a=la=l a 

We see that the gauge conditions JSJ, as well as (|24|l . are operator equations, which can be evaluated within 
commutators. We notice also that the definition i12-">[) of A is free of ordering problems, even though the 
commutators l|27l) are not canonical, but its components Aj in general do not close an angular momentum 
algebra, as shown by l|28b|l . 

In the quantum theory, a realization of the commutators (|27J) is obtained by defining P a as the projection 
of the gradient V Q on the gauge hyperplane 6 Q = 0, 



I Q 1 Id 

aJ a=l a f3=l 



These operators satisfy both (|27|l and the gauge condition l|24|l . They also satisfy relation (|25(l which, 
with p a j = —id/dr a j, is simply the chain rule for derivatives with respect to variables related by the 
transformation (|12|1 . Relations exactly analogous to (j27(l and (|29|l hold in Yang- Mills theories (compare (|29|l 
with the eq. between (6.13) and (6.14) in [H]). 



3.2 Quantum Hamiltonian 

The classical Hamiltonian in this gauge, (|26|l . was obtained from TCn in the gauge £ = by using the 
transformation l|25(l . The Hamiltonian operator can in principle be computed in a similar fashion, essentially 
by squaring l|25f) as an operator equation. That procedure works satisfactorily in the two-dimensional case 
PQ, but in three dimensions a more systematic approach is needed in order to handle the much larger amount 
of algebra required. The main difference between the two cases is that in two dimensions the operator A 
analogous to A in Ij25(l commutes with the Faddeev-Popov determinant [TJ , but that is not the case in three 
dimensions. Following we will first formulate the theory in terms of an appropriate set of independent 
generalized coordinates and their conjugate momenta. The results obtained in this intermediate step, which 
are of interest by themselves, will be transformed afterwards to the variables {R a } and {9 a }. 

The gauge conditions (JSJ) are defined by 9./V constants T aai , a = 1, 2, 3, a = 1, . . . , N, which constitute a 
set of three vectors F a with 3N components T aa i each, orthogonalized according to (|ll|l . We extend the set 
{F a }a=i to an orthogonal basis {To};^ of R 3N by arbitrarily choosing 3(N — 1) additional vectors {Tb}l= 4 
satisfying the orthogonality and completeness relations, 

N 3N 

^ m a T aaj T baj = y\ 2 a 5 ab , 1 < a, b < 3N, ° T aai r a j3j = m a 5 a p8 i: j , (30) 

q=1 a=l a 

which generalize pi)). We assume, for simplicity, that £H| — ... — = £H 2 , with $K 2 > an arbitrary 
constant. Furthermore, we choose all T aa i to have dimensions of length, so that S a , £2 a , and all have 
the dimensions of an inertia moment. We define a set of generalized coordinates q a , 1 < a < 3N, in the 
laboratory gauge £ = by, 

3N N 

r ai (t) = J2<la(t)^ i , q c (t) = ^c ai r ai (t) , 1 < c < W . (31) 

a— 1 ' a— 1 
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Similarly, we introduce 3iV — 3 independent generalized coordinates in the gauge (JSJ by. 

3N N 



Rai(t) = Qa(t)r aai , Qc(t) = ^£ T c m R m (t) , 4 < C < 3iV . (32) 



Due to the orthogonality relations (|30|l . the expression l|32|l for R a satisfies the gauge conditions JHJ- The 
dynamics in this gauge are completely specified by the 3iV independent variables {9 a }a=i an d {Qa} a =4i an d 
their conjugate momenta. Notice that the normalization of the coordinates q a and Q a is different. In (|32|l 
the Q a are chosen to be dimensionless, for later convenience, whereas in (I31|l the q a are defined so that the 
kinetic energy operator takes the simplest possible form, that of a Laplacian in Cartesian coordinates. 
The Hamiltonian TCn in the laboratory frame, Q, is given in terms of q a by 

1 3N d 2 

2 Q tt 9 ll 

from whence the expression for TIn in terms of {0 a } and {Q a } follows by means of a coordinate transfor- 
mation. The kinetic energy operator then takes the standard form of a Laplacian in curvilinear coordinates 
in configuration space. As shown in appendix 1X1 the result can be written as 

H = 1 V 9 « 9 

nN 2m. 2 J ^ dQ a dQ a 



E ^ + E w b A ») ^ m E^ + 1 ^w A/ 



(33) 



V 



2\A\J 

where |A| = det(A Qi ) with A ai defined in (JT3J and J = det(A/") 1/2 with 

^« = E i^c, . AA- 1 = ^ ^QT/QJ . (34) 

c=l c d=l 

The quantities £j a i with 4 < a < 3N appearing in (|33|) are defined as in eq. ©, of which they are an 
extension to a > 4. The inverse matrix 0" 1 , however, is defined only for a < 3. The expression <|33[1 for Ti^ 
depends explicitly on the constants T aa i with a > 4 through £2 a j with a > 4 and on the parametrization of 
U({6 a }) through AT" , both of which are largely arbitrary. Those dependences will disappear once we recast 
Hat in terms of P a and L. 

Using the relation, valid for any matrix depending on a parameter, 



6=1 

we obtain (compare (4.48) of [Sj) 

I <n A | A;,, 1 ) 

6=1 " " 6=1 



and hence, using 1|22|). 



1 1 6=1 6=1 



Therefore, we can rewrite TIn in terms of L as 

1 3N Ft Ft 
<H = 1 V 9 « 9 

nN VSPj^dQa dQ a 



a=4 

2 J U 2 ^ dQa m 1 ) ij 3 [ m 2 L> Qcj dQ 

\ a— 4 / \ c=4 



(38) 
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thus eliminating all explicit dependence of Hn on A^ 1 and the parametrization of U({9 a }). Applying the 
chain rule, from 1|32|) we have, 



d " d 



7^T = E r -7^' 4 < c <3^- (39) 



The fact that the R a are not independent, being related by (|HJ, does not affect l|39|l because of the orthog- 
onality relations (|30|l . Using the definition (|2fj|) of P a and the completeness relation l|30|) we get, 

3N 11 1 d N 

^-E^-— , ? ^ = £W«, 4<d<3A,, (40) 

and from l|4l7|l and completeness, 

1 3JV 1 » 1 3JV -. rj 

a=4 " c— 4 

with Aj given by (|25[) . Substituting 1)41) [) and (|4H into Hn in l|38|) we finally obtain 

N 1 1 
Hn = E o + - W^Jfa - A,) + V , (42) 

a—l 

in which all dependence on {r ace j}^ 4 has disappeared. The total Hamiltonian in this gauge is H = HN+H T t, 
with 7i rt from Q. Since s is a constant of the motion we can let X — > oo with s fixed, so that H T t vanishes. 
We are then left with an Af-body system described by Hn and the constraint L\ ) = —S\ ) = — U({8 a })s\ ) 
on the state space, with i = [Hn, s] = 0. We now turn to this constraint equation. 

3.3 Constraint and physical Hilbert space 

The wave function in this gauge ip({R a }, {9 a }, e) is required to satisfy the constraint J\if>) = (L+S)\ip) = 0, 
originating in the equation of motion l|4c|) . Expressing L in terms of pg a as in IL'l't , and using S = Us and 
(|13fl . we can write the constraint explicitly as 

i-j^- + \ ak Sk^ip{{Ra},{0a},e) = O , a =1,2,3. (43) 

We introduce the unitary operator U{{9 a }) 6], depending on {0 a } and acting on the Hilbert space of the 
rigid rotator, which satisfies the analog in that Hilbert space of H13|l . 

dlA 

—tf = iXakSk , a =1,2, 3. (44) 

Of a 

The matrix elements of U({9 a }) in the basis of eigenfuntions of s 2 , s z , (e\s,s z ) = Y 8S __{e), are the matrices 
D s s , s ^({6 a }) (given, e.g., in 20: in terms of Euler angles). Defining the physical wave function i/j({R a }, e) 

^{{R a }, {9 a }, e) = U({9 a })^({R a }, e) , (45) 

we see by direct substitution that (|45|l is a solution to the constraint equation (|43|) The wave functions 
ip({R a }, e) span the physical Hilbert space of the system. 

Some remarks about the form of the solution 145|) to the constraint are in order. By definition E = 
U{{0 a })e, and S = U({0 a })s. Since [si,(e)j] = iEijk(e)k, we have [Si,(E)j] = ie i3 k(E) k and from JTHJ), 
[Lfc, (E)i] = —i£kin{E) n . From these relations, and taking into account that s is a differential operator 
on the Hilbert space of the rigid rotator, we obtain U({8 a })x(e) = x{E) for any x(^). Thus, (|45|l can be 
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rewritten as ip({R a }, {9 a }, e) = ^({i? Q }, E). Had we formulated the theory in terms of {R a }, {9 a } and E, 
the momenta conjugate to 9 a would have been pg a = —A a iJi, instead of (|21|) . and the constraint Jtp = 
would have led to ip n °t depending on {9 a }, ip({R a }, {9 a }, E) — ip({R a }, E), the same result as 1)450 . 

From (|44l) we obtain [Lk,U] = —SkU : which leads to U^LkU = L k —U^SkU- Either from this last 
equation or from the constraint, on the physical wave functions we have IvLyUtj) = —Iv SkUip = —Skip- 
Therefore, in the physical Hilbert space the Hamiltonian (|42|) acquires the form, 

N 1 1 
H N = U^H N U = o 7f P <xiJP ai + W7f(si + A l )M^ 1 J{s J + A,-) + V . (46) 

Similarly, we can write TIn in terms of the independent coordinates {Q a } and their conjugate momenta, by 
just substituting — s for L in i)38|) . Thus, once the constraint has been solved the angular variables enter 
the dynamics only through the dependence ofH.pf on s, the angular momentum of the rigid rotator in the 
lab frame. Notice that although s is a constant of motion in the full Hilbert space of the theory, in general 
[s,ii.pf] 7^ because [s,U] ^ 0. Thus, although s 2 can always be diagonalized simultaneously with Ti^ 
within the physical subspace, s z in general cannot be diagonalized. The physical reason is that within the 
physical subspace the matrix elements of s are equal to those of — L, which is not conserved. 



3.4 Inner product in Hilbert space 

In order to find the inner product in the gauge 6 a — we transform its expression (JjJ in the gauge £ = 
by means of the Faddeev-Popov technique \F3\ . For that purpose we first find an appropriate resolution of 
the identity over the group 5*0(3), fixing on the way any Gribov ambiguities ^212] inherent in the gauge 
conditions. The invariant integration over 50(3) is given by (see, e.g., |21|) 

3 

/ dgf(g)= ff[d9 a \A\f(g({9 a })) , (47) 

where on the l.h.s. the integration variable g takes values in 5*0(3) and / : 50(3) — > C. On the r.h.s. of (|47|l 
the integration extends over all of parameter space and |A| = det(A a j) with A a i defined in 1)13)) . From (|36|l 
and gZI), the operators L defined in ^ are hermitian, J lf a=1 d9 a \A\ip* L k <p = {J l\ 3 a=1 d9 a \A\cj)* L k ip)* . 

3.4.1 Resolution of the identity. Singularities of the coordinate frame (Gribov ambiguities) 

With the integration measure 1)47)1 . the resolution of the identity for the gauge 1)10)1 takes the form, 

1= ff[d9 a \A\f[S l^& b ({U(9)r a })) J ({U(9)r a }) 9 (J {{U(9)r a })) Q(f) , (48) 

J a=l 6=1 VJth J 

where the gauge conditions are conventionally written as & b /9K b , and J({r a }) is defined after (|3*3"|) . To 
obtain |@SJ), let 9 = {9 0a } be a root to & b ({U(9 )r a }) = 0, b = 1, 2, 3, for some fixed {r a }. Then, 



JY 



& b ({U(9 + S9)r a }) = m a T baj 5r a j , 6r a j = y^ y 59 c -^-(9 )r ak = ^S9 c A cn Sj n iUi m (6o)r am (49) 

a— 1 c— 1 ' c— 1 

where for the last equality we used 1)13)) . Using ©, we rewrite 1)49)) as 

1 3 1 

— e b ({U(9 + 59)r a }) = S9 c —A cn Q bn ({U(9 )r a }) , (50) 

c—1 

so that, 

det ( *&£**L ) = |A| det ( ^-Q bn ({U(9 )r a })) = |A| det (Q bn {{U(9 a )r a })) , (51) 



ii 



and therefore 

ft * ({U(0o)r a })) = £ ft W» - **) , (52) 

where the sum extends over all roots 9o. In order for the l.h.s. of (|48|l to be 1, 9o must be unique. The 
linear gauge conditions, however, have in general Gribov ambiguities leading to a discrete set of roots 9q. We 
assume that we have chosen the parametrization U({9 a }) so that |A| > 0. But det(£2) can vanish at some 
configurations {r a } at which the gauge conditions are singular. Given a configuration {R! a = U({9' a })r ai 9' a } 
with det(O) < 0, we can always find a gauge-equivalent one {R a = U({9 a })r a , 8 a } with det(£2) > 0. Thus, 
we restrict ourselves to those configurations satisfying, 

< = (det(AO) 1 / 2 = J . (53) 

If this supplementary condition were enough to remove all ambiguities, together with (|52"|l it would lead to 
1148( 1 . Unlike the two-dimensional case pQ, however, choosing the sign of the Faddeev-Popov determinant J is 
in general not enough to remove the ambiguities. Further supplementary conditions may be required which 
are of the form Fx > 0, . . . , F r > 0, where the Fj are r functions of the particle coordinates (as many as 
necessary to fix the gauge), such that for each j it is true that Fj = implies J = 0. This set of additional 
supplementary conditions is symbolized by the factor 0(.F) in l)48[l. 

A simple example will illustrate the previous discussion. Assume that for a system of N > 3 particles 
we want to choose a coordinate frame rotating so that particle 1 is on the X axis and particle 2 is on the 
X—Z plane for all t. That frame is not well defined if particle 1 is at the origin or particle 2 is on the X 
axis. The gauge conditions defining the frame are Si = Rxy = 0, &2 = Riz = 0, 63 = R2Y = 0, leading to 
dct(£2) = —R\ x Riz- As expected, det(O) = if Rxx = (1 is at the origin) or R 2 z = (2 is on the X 
axis) . Those singularities stem from the fact that there are four ways to choose the rotating frame, depending 
on whether we choose R\x < 0, Riz < for all t. By requiring J > we must have Riz < 0, which fixes the 
ambiguity only partially. In order to completely fix the gauge we have to impose a supplementary condition 
such as Fx = Rxx > 0. Clearly, Fx = implies det(£2) = = J . Alternatively, we may exploit the fact 
that if we do not impose the condition Rxx > every configuration is counted twice (except for those with 
Rix = which are counted once, but they have zero measure and do not contribute to l|48|l). Thus, in this 
example we may omit the factor Q(J-) on the r.h.s. of (|48|l and set the l.h.s. to 2. 

The general case is analogous to the simple example above. We either have to include in (|48l) the factor 
appropriate to the gauge conditions, or replace it by a factor 1/(1 + N({R a })) with N({R a }) the 
number of gauge-equivalent copies of each configuration {R a } satisfying the gauge conditions and J > 
|llj . In those cases in which N({R a }) is a constant over all of configuration space (except maybe for a 
zero-measure set) we can omit those factors, absorbing them in the normalization of the inner product. 

3.4.2 Inner product in Hilbert space 

The inner product in this gauge is straightforward to obtain from J7J by using the Faddeev-Popov trick with 
the resolution of the identity (|48|l . We briefly sketch the derivation in order to highlight the relationship 
among wave functions in the gauge £ = and those in this gauge. We rewrite Q as, 

N 

m) = K / J] d 3 r' fj d 2 e (4>M({U(e a )r' }, e) , (54) 

J 3=1 

where k > is a normalization constant to be chosen later (in JJJ, k = 1), and we performed a change of 
variables r'p = U(9 a )rg with U(9 a ) an orthogonal matrix. (In 154(1 we temporarily introduced the notation 
"00 for wave functions in the gauge £ = for convenience.) Inserting l(48|) in 1(541) . exchanging the order of 
integration, and changing variables back to rp, we get 

(M) =« / II^I A I / IPV^IP (w &b ) J®(JMF)&M({r },£) , (55) 

^ o=l /3=1 6=1 V 6 / 
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where we omitted the argument {rp} in <5 and J for brevity. We can now choose n to be the reciprocal of 
the volume of the rotation group. Identifying, up to a phase factor, the physical wave function -0({_R a }, e) — 
ipo{{Ra}, e) for {R a } satisfying the gauge condition we finally obtain, 

N 

(<t>W = I TT d 4 Rpd 2 e \\S[ --e b ) je(J)e(F)(<p*i>)({Rf3}, e) . (56) 



0= 

This gives the inner product in the physical Hilbert space of the system. The matrix elements for the 
Hamiltonian can be computed with the operator Ti.^ of l|46|) and the inner product (|56|l . The hermiticity of 
7ijv with respect to (|56H follows by partial integration, taking into account that P a i and Aj are homogeneous 
first-order differential operators, with constant coefficients, satisfying [P a i, &b] = = [Aj, @b], and that 
J8(Fj) = since, as mentioned above, Fj = implies J = 0. 

We remark also that for states satisfying the constraint [L + S)\ip) — 0, using the definition (|45|l of the 
physical wave function and W SU = swe get the equality 

J J{de a \h\ [dii<fr%{R a },{9 a }),tyL i il>{{R a },{0 a }),ty = ± J d^r({R a },e)(~ Si yj({R a },e) , (57) 

a—l 

where we denoted by d\x the measure appearing in l|5t)l) . The factor 1/k appears in (|57|) due to the normal- 
ization we chose for the inner product in the physical subspace. From (|57|l we see that for physical wave 
functions the operator —s gives the matrix elements of L. 



3.5 Reduced Hamiltonian and Weyl ordering. Quantum potential. 

When working in curvilinear coordinates it is often convenient to redefine the state space by absorbing the 
Jacobian in the wave functions, thus eliminating it from the integration measure in the inner product and 
from the kinetic energy operator. That is the case, for instance, when a perturbative expansion of J contains 
terms of many different orders. Furthermore, the reduced Hamiltonian is easier to cast into Weyl-ordered 
form, in which the relation between the operator and path-integral approaches is straightforward. 

The Hamiltonian l|4t)|) has been simplified by restricting it to the physical Hilbert space of gauge- invariant 
wave functions "0, so it is not of the form of a Laplacian in curvilinear coordinates. Thus, we go back to the 
form Ij33|l for T~Ln, in which the angles {9 a } and their conjugate momenta pea = —id/d6 a appear explicitly. 
From i|33|) and (|A.17|) we get, 

i ^ d d l ^ a ^ , d d \ 

- ^gg fa A »^>ik*k) w - { A ^ QbkM *wA) w (58) 

where Vq is the quantum potential given below (see Q65[)) and (• • • )w indicates Weyl-ordering (e.g., IjATlj) - ). 
Let us introduce the notation 

3JV 

^ij = £ ijkRak = ^ ^^aaiQaj , (59) 
a—l a 

where the second equality follows from and 1|30|) . With this definition and relations l|40fl and (|24ll . we 
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can rewrite ()58|l in terms of P a operators as 

N „ , N 3 



1 1 1 / 

a = l Q,/3=l 6=1 V 

a— 1 a.b— 1 x ,! 



(60) 



This operator is to be applied to wave functions of the form ip({R a }, {6 a }, e) = \A\ 1 ^ 2 J' 1 / 2 'ifj({R a }, {0 a }, e) 
or, using (|4*5l) . 



^({-R ct },{^} ) e) = |A|=W({^})^({H a },e) with e) = ^({i^}, g) . ( 6 l) 

i/j is the reduced form of the physical wave function tp of 145I) . 

From H36f) we obtain the form of the angular momentum operator on the first line of (|60ll when applied 
to wave functions of the form (|61fl 

S 5 ( A ^7^ + \w b A ^) * = -l A l^^({^i)S = \kM{e a }) Sk 4 , (62) 

6—1 

the second equality following from the discussion immediately above (|46|l . After appropriately rearranging 
the angular operator on the second line of l|6(J|) we can apply (|62[1 to it as well and, taking into account that 
N^ 1 is symmetric, we get 

a, 6=1 



de b de a v Ja de a de a > a j y kb de b de b 

d0 b 86 a V T 



a,b=l 



Thus, gathering fi ty . and (E3 together we obtain 
f W 1 1 N 



a,/3=l 



2 ^— ^ V w 1x1 jk w j 2 J 8 Jfc ^ <90 b <90 a 

a=l a, 6=1 

(64) 

The quantum potential is computed in appendix IA. 21 with the result, 

1 3 <9A _1 <9A _1 1 N 3 

V Q = E + Vl + V 2 > V i = -gEE ^Qi,lT cal £l^ d T dak eu- P e P i m (65) 

a, 6=1 a=lc,d=l 

- 1 N -1 ( 3 i \ ( 3 i \ 

(8,7=1 V 0=1 /V 6=1 / 



14 



(Notice that both Vi,2 are 0(|il a | 2 ) as \R a \ — > oo). Thus, the Weyl-ordered, reduced Hamiltonian H jv 
= J 1/2 H N J~ 1/2 = WH N U acting on the space of reduced physical wave functions ip ({R a }, e) is given by, 

= E 2^ P " + ^ ^ (E P /^« + ^ (E + ) ) + Vi + V 2 + V , (66) 

where the second term is, explicitly, 

N 




' 3,7=1 

1 N 1 



(67) 



4 

0=1 



Given two states |</>) and |^>) represented by the reduced, physical wave functions <f>, ip, their inner product 
is, according to (|5l)j) and (51), 



TV 3 

MV>> = / fl rf 3 ^rf 2 en^f^6 b ) e( l 7)6(^)0*({^},e)^({^},e) . (68) 

J 13=1 6=1 V " nb / 

The matrix elements of TL n computed with this inner product are, of course, identical to those of the 
operator Hjv of (|4*fi)l computed with the inner product (55) . 

As pointed out before, there is a close formal analogy between the results given above for the quantum 
theory in reference frames defined by linear conditions and the corresponding results in Yang-Mills theories 
in non-covariant linear gauges. The derivation of the Hamiltonian given here parallels that of [Sj- Thus, 
the quantum potentials Vi,2 from l|65|) are formally analogous to the corresponding expressions (6.12) and 
(6.14) in |5]. The kinetic energy operator in H46fl and in the Weyl-ordered form H66|) are formally equivalent 
to (4.62) and (6.15) of J5|, respectively. In order to make the formal analogy clear, we notice that the space 
derivatives appearing in field-theoretic expressions must be mapped to zero in the mechanical case considered 
here. Thus, Q a i from © is the analog of the expression T^Dk in the notation of [5j- From (|25|1 and l|59|) . 
A; = — Yj a =i m a^ > ijPaji which is the analog of S^P/ ~ —P-®i and, similarly, Af^ 1 from is identified 
with (T^r^TjT^lTl)-' in 0. 



4 Center of mass motion 

In this section and the next one we set U = in the Lagrangian and take into account the translation 
invariance of Q in order to separate the center-of-mass degrees of freedom. Since the center of mass motion 
is dynamically trivial, we restrict our treatment to dynamical states with vanishing total momentum. 
The Lagrangian Q is invariant under time-independent transformations of the Euclidean group, 

r' a = Ur a + u , e = Ue , (69) 

with U an orthogonal matrix. We define the covariant derivatives 

D t r a =r a -£r a - p , D t e = e - £e . (70) 

Under time-dependent transformations r a and e transform as in (53) and, 

£' = U£tf + Utf, p' = Up + u-$'u, (D t r a )' = UD t r a , (D t e)' = UD t e. (71) 
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Substituting time derivatives by covariant ones in we obtain a Lagrangian invariant under the time- 
dependent transformations (|69fl and l|71|l . 



^ N N 



a=l 



where £jv and £ r t have the same form as in 0. The equations of motion for r a and e take the form (0} 
when expressed in terms of covariant derivatives, but in this case the derivation is slightly more involved 
because now D t r a as given by (|7U)l does not depend linearly on r a but, rather, affinely (see appendix |B|) . 
The angular momenta I and s are still given by (JSJ, with D t r a from J7DJ. We define Z cm = Mr cm A D t r cm , 
where r cm is the center-of-mass position, D t r cnl — r cm — £r cm — p, and M = X^=i TO a- As shown in 
appendix [5] the equations of motion lead to, 

" t) (I - i cm ) = = ( ± ~ t) s , (73) 



dt v v ^' \dt 

which is the same as (|SJ), with (I — Z cm ) instead of I. Thus, the magnitudes of (Z — Z cm ) and s are both 
conserved and frame- independent. Due to the fact that we are now including (time-dependent) translations 
as symmetries of the system, the role played by I in sections [5] and [3] is now played by (I — l cm ). The eqs. of 
motion for £ and p are now of the form I + s = and p cm = (with p cm = MD t r cm ). Thus, in particular 
^cm = and, like in sectional the total angular momentum of the system vanishes. 

If we choose the gauge conditions £ = = p, corresponding to the laboratory frame, we recover the 
Lagrangian l|T[l. constrained by the eqs. of motion for £ and p in this gauge, 

N N 

m a r a A r a + Xe A e = I + s — , m a r a — . (74) 

a—l a—1 

These constraints are first class. In the quantum theory they restrict the state space, (I + s)ip = 0, 
J2 a =i ^ ai> — 0' analogously to Gauss law in Yang-Mills theories [§]. Except for the additional constraint 
on the center-of-mass momentum, the quantization in this gauge is carried out exactly as in section \2. II 
We can now proceed along the same lines as in section imposing on the system the gauge conditions 



I N 

& a ({R a }) = 0, a = l,2,3, C({R a }) = — J2 m P R f3 = °, ( 75 ) 



0=1 



with & a defined in (JHJl. I|75JI defines a reference frame in a particular state of rotation, with origin at the 
center of mass. Like in scction[31 in the rest of this section we denote vectors referred to this frame by capital 
letters, while lower-case symbols denote lab frame quantities. The gauge conditions l)75|) are not mutually 
consistent unless & a are translation invariant, 



N 



2_^m a T aaj =0, a = 1,2, 3, (76) 



the condition C = being clearly rotationally invariant. Furthermore, we assume that 6 a satisfy (|10|l 
and (H) . The gauge transformation from the gauge £ = = p is of the form (|69|I -H71 |I . with parameter 
u = —Ur cra , where r cm = $3„_i Ti a /Afr a is the center of mass in the lab frame, 

R a = U (r a - r cm ) , E = Ue , £ = UU f , p = -Uf cm . (77) 

The transformation 17711 mixes the particle degrees of freedom R a with those of the center of mass and the 
rigid rotator, just like (|12J) did in the non-translation-invariant case. In particular, in these variables pg a is 
not linearly related to L. That mixing is avoided, as in section |21 by trading the dynamical variables 
E, £, p for {R a }, e, {8 a }, r cm . Substituting the last three of (|77|l into the Lagrangian (f72l we get, 

C=C N + C lt + £cm , £cm = , (78) 
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with Cn and £ rt now given by l(T7|) . The Lagrangian (|7S|) is supplemented by the gauge conditions (j75|l 
holding as strong (operator) equations, and the constraints J = L + S = and p cm = Mr cm ~ valid as 
weak (state space) equalities. 

We keep the definitions (|13[1 of A a i and \ a i from section [3J The expression (I14f) of £ in terms of {9 a } then 
holds unchanged since £ = UU^ just like in sectional From C in (|78|l we can then derive the relations (|21|l 
for pg a and the angular momentum L in this gauge. The classical expression (12111 for L, in turn, together 
with the transformation law (|77J) lead to the relation 

L = U(l- l cm ) , (79) 

where the center-of-mass angular momentum in the lab frame is defined as Z cm = r cm A p cm = Mr cm A r cm . 
From the relation l|77[) between R a and r a we can derive an expression for dU / dr a jU^ by following the 
same steps leading to (I16f) in sectional The result is that (I16|l remains valid without modifications and that, 
due to the translation invariance condition (|76() for & a , U does not depend on r cm , 



a—l 

In particular, [Z cm , U] = 0. From Ijltifl . in turn, the commutator (|18fl of I with U follows. Thus, with the 
commutator {THJ, the transformation relations l(77|) and {73, and (JSUJ), we recover all of the commutators 
(I19fl and also, 

Furthermore, [Z C mi, s^] = 0. Thus, since J = L+S = U(l— i cm +s), with 1|81|) we find [Z C mi, Sj] = = [Z C mi, 
and then all of the commutators (|20|) follow. In summary, with the exception of eq. (|79|l . all of the results 
of section \'S . II remain valid in this case. 

The relation among linear momenta in the gauge £ = = p and the gauge 6 a = = C analogous to 
(|25[l takes the form, 

/ 3 \ m 

Pcy = U kj iP ak + '^2m a T aak £l~a(L n - A„)J + -j^p cmj (82) 

with A defined as in (|25|l . Correspondingly, the classical Hamiltonian is given by (|26|l with the addition of 
the center-of-mass kinetic energy p^ m /(2M). Due to the additional gauge conditions <L = the fundamental 
commutators Q27JI become, 

[R ai ,P/3j] = i ^SapSij - -j^Sij - X] ^§- r <"Mra/3j^ , (83) 

and [P a i, & a ({Rp})] — = [P a i, Cd-R^})]. The differential operators realizing this algebra are obtained in 
the same way as those in (|29H . which is now modified to 

~ * dR ^ h m * am h 1 dR ^ M h 1 dRpi ' 

These operators satisfy & a {{Pa/m a }) = = J2a=i ^ a They also satisfy JB2Jl, with p a j — l/id/dr a j. The 
additional term in (|84|l with respect to l|29|l does not modify the form of A as a differential operator. Clearly, 
the commutators of R a and P a with r cm and p cm = l/id/dr cm vanish strongly. 

4.1 Quantum Hamiltonian 

The Hamiltonian operator is obtained by the same procedure as in section T3. 21 with obvious modifications. 
We omit all calculational details and quote the results only, after establishing the appropriate notation. 
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Like in section we extend the gauge coefficients T aa i, a — 1,2, 3, a = 1, . . . ,N, to an orthogonal 
basis of M. 3N . That is, we consider an extended set of coefficients T aa i with 1 < a < 3N, a = 1, ...,N, 
i = 1, . . . , 3, satisfying the orthogonality and completeness relations (I3U|I . For simplicity, like in section HOI 
we set the normalization constants in (|30|l to be independent of a for 4 < a < 3N, 9l 2 = 9l 2 > with 
9i 2 a constant at our disposal. Furthermore, for 37V — 2 < a < 3N we choose the coefficients T aa j to be 
independent of a. Specifically, we set, 

T aal = —S (a _ m+3)i , 3N-2<a<3N , l<a<N . (85) 
V M 

With this choice the orthogonality relations (|30|l with 3N — 2 < a < 3N and 1 < b < 3N — 3 read, 

N 

m ^baj = , 1 < 6 < 3N - 3 . (86) 

a=l 

In particular, (|86|l contains the translation invariance conditions (|76|l for (Sf,. 

We can now define the generalized coordinates {q a }a=i by (|31|l . We see that for 3A^ — 2 < a < 3iV, 
q a = V r /V/r cm ( Q _3jv+3), i-C, up to a normalization constant the last three q a are the components of r cm . The 
analog of 1-iLl is now, 

3JV-3 N 
Rai(t) = Qa(tWaai , Qc{t) = ^ ^^c m R m (t), 4 < C < 3iV - 3 . (87) 

a=4 a=l 

From H30(l and Ij85(l the expression (|87|l for R a satisfies the gauge conditions (|75|l . The dynamics in this 
gauge is then completely specified by the 3N independent variables {0 a }a=ii {Qa}a=4 3 and r cm , and their 
conjugate momenta. In those variables TLn is given by (|33|l with only two modifications: first, the sums over 
indices running up to 3iV now run only up to 3./V— 3, and second, the addition of the term — l/(2Af)<9 2 /9r cm 2 . 
The definitions (|34|l of Afij and its inverse, and of J and |A| remain unchanged. Similarly, TCn is expressed 
in terms of R a , their conjugate momenta P a , and L by (|42(l . but now with the momentum operators P a 
from (|84JI , and with the addition of the center-of-mass kinetic energy term. 



4.2 Physical Hilbert space. Inner product. Weyl-ordered Hamiltonian 

The wave function in this gauge ^({H Q },{0 a },e,r cm ) is required to satisfy the constraints p cm ijj = and 
(L+S)ip — originating in the equations of motion for p and £ from the Lagrangian l|72() . The first constraint 
is trivial to solve. Considering wave functions ip independent of r cm , we are left with the constraint on the 
angular variables which, since [p cm , L] = = [p cm , Sf], can now be treated exactly as in section 13.31 Using 
the same notation as in (|45|l . the solution to the constraint equations is of the form, 

i/>({R a }, {9 a }, e, r cm ) = U({0 a })4>({R a }, e) . (88) 

Within the subspace of physical wave functions ip({R a }, e) the Hamiltonian Hn = U^TInU is given by l|4^|) . 
with the momentum operators P a from (|84|l . 

The discussion of the inner product from section 13 . 41 rea uires only minor changes in order to adapt it to 
the translation-invariant case. Besides the resolution of the identity (|48|l for the rotational gauge conditions 
& a = 0, we have to fix the translational gauge by means of a resolution of the form, 

3 

1= / ^uHS^dr^ + u,) . (89) 

i—l 

Inserting this factor of one together with l|48|l into the canonical inner product {7|), we obtain (4>\ip) in terms 
of wave functions in this gauge. A technical detail is that, after applying the Faddeev-Popov procedure, 
the volume of the symmetry group appears as a prefactor in (4>\ip) (see <|55|l ^) . In this case, the volume 
of the translation group is infinite so an appropriate limiting or regularization procedure must be applied. 
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Assuming that has been done, the resulting inner product in terms of physical wave functions analogous to 
is 

W#> = fli d 3 R a d 2 ef[s(^& a ) 5® (£)e(F)Q(J)J($*$)({R a },e) . (90) 

a— 1 a— 1 ^ ' 

The factor 0(.F) is exactly as discussed in section l3.4.1l since the center-of-mass condition does not introduce 
further Gribov ambiguities. 

We consider, finally, the form of the Weyl-ordered reduced Hamiltonian. Taking proper account of trans- 
lation invariance as described above, the analysis of section 13.51 remains valid mutatis-mutandis. Defining 
the reduced physical wave functions 

^({R a },e) = J^({R a },e), (91) 

with as defined in Qgft, the reduced Hamiltonian H N = J X I 2 H N J- X I 2 = J X I 2 WU N UJ- X I 2 is given 
in Weyl-ordered form by (|66|l . with P a given by l|84|l . The quantum potentials Vi,2j i n particular, are still 
defined as in JB5J). The reduced inner product is immediately obtained from lj9"5jl and (foTT) . 



5 Quasi-rigid systems in the Eckart frame 

We assume now that the potential energy V (with U = 0) has a minimum for some configuration {z a } of 
the system, such that Vo = V({z a }) < V({r 7 }) for all configurations {r 7 }. Due to the invariance of V under 
the Euclidean group E$ any configuration {z' a } related to {z a } by a transformation of the form i|69[l is also 
a minimum. Denoting by M\> the manifold of configuration space defined by V({r" 7 }) = Vo, we assume that 
the quotient My/E^ is a discrete set. The configurations of minimal potential energy are therefore rigid. In 
this section we discuss the quantization of the small oscillations of the system about these rigid equilibrium 
configurations. We will denote by {Z a } the unique (up to discrete degeneracy) minimum of V satisfying, 

N N 

^m a Z ai Z aj = , i^j, ^m Q Z Q =0. (92) 

a— 1 a— 1 

We will restrict ourselves to considering only systems for which the inertia tensor for the equilibrium config- 
uration {Z a } is non-singular. The small oscillations of the system are described by trajectories of the form, 

r a {t) = z a (t) + dr a (t) with z a (t) = U{t)Z a + u (93) 

for some orthogonal matrix U(t) and u appropriately chosen so that 8r a (t) are small with respect to their 
characteristic scale for all t. Since we restrict ourselves to states with vanishing total momentum, the 
translation vector u in l|93|> must be time-independent. 

It is convenient to apply the inverse of the gauge transformation defined by the second equation in Q93Jl 
in order to switch to a reference frame, the "body frame" of the rigid equilibrium configuration, so that 

r a (t) = Z a + 6r a {t) . (94) 

This fixes the gauge only to leading order in Sr a . We fix the residual gauge freedom by imposing a gauge 
condition on Sr a , which amounts to correcting the definition (|92() - (|94|l of the reference frame by small 
quantities of first order. We choose the origin of the reference frame at the center of mass, so to first order 
in Sr a the gauge conditions must be of the form (|75|l . The choice of the coefficients T aai is arbitrary as long 
as (I76|l is satisfied. We then have, 

R a (t) = Z a + SR a (t) , & a ({SR a }) = , €{{6R a }) = . (95) 

The Eckart frame corresponds to choosing T aai = e a jiZ a j, a = 1, 2, 3 |T21I9]. With this choice the normal- 
ization constants 9\ 2 , a = 1,2,3, of (|11|) (and also (|30|l and (JEUO are given by the inertia moments of the 
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equilibrium configuration, and the matrix Q({R a i}) of JjJJl i s J 

N 

Q ai ({SR a }) = m 2 a 6 ai + 5Q ai ({8R a }) , 5Q ai ({8R a }) = ]T m 1 {Z 1 ■ dR^5 a i — SR^aZ^i) , a = 1, 2, 3 . 

7=1 

(96) 

Thus, since JH^ 7^ by assumption, for small 6R a the condition det(0 ai ({i? Q })) ^ is satisfied. 

The momentum operators P a and their fundamental commutators are as given in i|83[l and 18411 . with 
the replacement of d/dR a i by d/d8R a i. The operator A defined in i|25|) can be rewritten in the form, 

a, = £ e , fc «? QJ i^ - E ^ E ^ feZ -7^b ■ (97) 

a=l a=l " a=l 

As expected in this gauge |^|Q], its coefficients are of 0(SR a ). Similarly, the operators Yl a =i ^ij^ai and 

Y2a=i PaiDfj appearing in H n are of 0(SR a ). The Hamiltonian Hn, given by Ij66|) with P a from l|84|) . is 
obtained perturbatively by expanding Ij66|) in powers of 5R a . From this point of view, the elimination of 
the Jacobian J from the kinetic energy, as indicated in sections 13.51 and 14.21 is particularly convenient in 
perturbation theory. The inner product, finally, is given by ()90|l with the modification (|91|) . and with SR a as 
integration variable. Once the equilibrium configuration {Z a } has been chosen, its body frame is uniquely 
fixed. The Eckart frame is then equally well defined as long as {SR a } are small. Thus, except in those cases 
in which the equilibrium configuration {Z a } is exceptionally close to a zero of J ', we can neglect the factor 
Q{J) in (|90[) since large displacements 8R a which could drive J to zero should be exponentially suppressed 
by the wave function. Similarly, we also expect to be able to neglect in (|90[l in perturbation theory. 

As a minimal illustration and consistency check of the formalism we analyze below a simple example with 
N = 3, and briefly comment on the N — 4 case. The natural variables for quasi-rigid systems are normal 
coordinates, so below we recast the Hamiltonian in terms of those coordinates. We remark, however, that 
the results of the previous sections are more general than the simple examples considered here, and can be 
applied to non-quasi-rigid A-body sytems, both in the operator and path integral formalims. 



5.1 A simple example with N = 3 

The simplest possible model, within our assumptions, consists of three particles of equal mass m interacting 
through a two-body potential V as in with U = and V a p = V independent of a, (3. V{r) is assumed to 
have an absolute minimum at r = a > 0. The classical equilibrium configurations arc then those in which the 
particles lie at relative rest on the vertices of an equilateral triangle of side a. An equilibrium configuration 
satisfying (|92|l . unique up to permutations of the particles and discrete rotations of the coordinate axes, is 

z '="H'-57r°) ■ ^"G'-stt ) • ^"("'Tr") ■ (98) 

with the inertia tensor ma 2 /2 diag(f , 1 , 2). 

The Eckart gauge is defined by (|75|l with T aa i = e a jiZ a j, a = f,2,3, which are normalized to = 
ma 2 /2 = 9^, d\\ = ma 2 . Those gauge conditions make the planar nature of the problem apparent, since 
they imply that 5R a 3 = = P a 3, a = f,2,3, as operators. This leads, in particular, to the operator A of 
(|25|l having two null components Ai = = A2 (as operators), with A3 being conserved and having integer 
eigenvalues, as shown below. Thus, in this example A is an angular momentum operator, though two- rather 
than three-dimensional, and can be rightfully termed "residual" angular momentum as in pp. The three- 
dimensional rotations of the system are taken into account by the total angular momentum operator s. 

Setting V" (a) = muj 2 , the quadratic terms in an expansion of V about {Z a } are, 

= ^ E (-{Z a ~Z p )-{8R a -8R [3 )\ . (99) 

Using the gauge conditions we could eliminate six degrees of freedom, describing the system in terms of, e.g., 
SRix, 8R2X, 3R2Y, and their conjugate momenta. A better approach is to use a set of normal coordinates 
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{5Q a }a=4 as discussed in sections [3.21 and |4~T1 Thus, with T aai , a — 1,2,3, as defined above and r aQi , 
a = 7, 8, 9, as defined by (|85|l . we can choose T aa i, a = 4, 5, 6, to be those eigenvectors of the quadratic form 
associated with V(2) which satisfy the orthogonality conditions (|30|l (in particular, (|86|l '). and normalized 5 
to £R 2 = h/uj. Those T aa i, a = 4,5,6, are the vibrational normal modes of V( 2 ), whose associated normal 
coordinates 8Q a are given by JS7J, 



Similarly, from H39|) . I|4UI) (with 3A*" — 3 instead of 37V) we obtain the relation between d/dSQ a and either 
d/dSR ai or P ai . The result is given by (|100fl with <5Q a substituted by l/id/dSQ a on the l.h.s., and SR a i 
substituted by either h/(muj)l/id/dR a i or l/(mu>)P a i, respectively, on the r.h.s. 

The quadratic piece of the Hamiltonian Ti n (henceforth Ti ) of section l4~2l can be written as, 

8.-^£l+>to-T£(-rag I +'**). <^-|— a. -3-«- (">u 

/3=1 a=4 \ / 



7io is the lowest-order Hamiltonian in a perturbative expansion in powers of e = ^Jh/(muja 2 ) <C 1. From 
the definition (|25|l . or equivalently from H97[). we find the residual angular momentum as, 

:iA " 3 SQai Id , , , „ , 1 / a _ 9 



A '=E^-7iI7r and therefore = = A 2 , A 3 = - [SQ 5 —- - SQ 4 —- , (102) 
^ i ddQa i \ oSQi oSQ 5 J 

with <50 M defined in Qjfft. To C(e 2 ) the quantities entering ft are ml = SH| = 9^§/2 = /i/(2cje 2 ), Vv^ 1 = 
l/9^ 2 -- ) (5( i )j + 0(e 3 ), Vi = C(e 3 ), and the anharmonic terms in V, which are 0(e 3 ). The expansion of 
\>2 = est. + 0(e 3 ) starts at C(e 2 ) but the lowest-order term is a constant, which we drop. Retaining only 
terms through C(e 2 ) in Ti , and expressing them in terms of normal coordinates, we obtain 

we 2 ( o o 1 



H=H + Hx, Hi = ^- ^ 2 -^ + ^( S3 + A 3 ) 2 J + 0(e 3 ) . (103) 

where we used A3 as defined in Q102|l . and dropped all constant terms. The operators s 2 , S3 and A3 all 
commute with each other and with the Hamiltonian. The physical meaning of the Hamiltonian \ 103(1 is 
apparent from l|57|) : for physical wave functions the operator s gives the matrix elements of —L so that, to 
this order, (|103fl corresponds to Ti\ — l/2AfH (Li — Ai)(Lj — Aj), with M^ 1 the inverse of the equilibrium 
inertia tensor and L A the difference of the total and residual angular momenta in the Eckart frame. 

Since we are not going to compute beyond 0(e 2 ) in perturbation theory, it is easier (to this order) to 
solve the eigenvalue problem for (|103fl exactly rather than as a perturbation about the Hamiltonian (|101|) . 
We introduce cylindrical coordinates in the space of SQ^ ^^ 

5Q 4 
SQ 5 



P = JSQj +SQI, <p= arctan ( ^ ) , C = SQe , (104) 



and classify the Hamiltonian eigenfunctions and eigenvalues according to the eigenvalues oiTio, A 3 , s 2 and 
s 3 (with quantum numbers denoted by (n,n^), A, t, m, resp.). The wave functions are, 



"" (p,<p,C,e) - R nW { P )$> x {v)Z n( {QY im {e) 




R nW { P ) = pWLW ( \/V ] e-vW/2 ^ Q x{ip) = t ZnAC) = Hn({V z 0e -3C/2 



(105) 



3 In this section and the following we restore h in all expressions. 
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where we omitted a normalization constant, n, tiq and £ are non-negative integers, A and m are integers, and 
the spherical harmonics Yg m , associated Laguerre polynomials L\ and Hermite polynomials H n are defined 
in the standard way in quantum mechanics 20 . The dependence of the wave function l|105|) on e only 
carries the representation of the rotation group appropriate to a state of angular momentum £. We could 
as well suppress the dependence on e and define the wave function to be a column with 2£ + 1 components, 
depending only on the three vibrational variables p, (p, £. The energy eigenvalues are, 

rplm _ p (0) , 

E ( °L ( = (ji c + I) + Vf( 2 " + l A l + !)) > 41 = ^ 2 (V + 1) - ™ 2 + \(™ + A) ; 

(106) 

In the vibrational ground state n = = A = 0, i|106fl reduces to the spectrum of an axis-symmetric top. 
5.2 A simple example with N = 4 

Adding one more particle of the same mass to the model of the previous section we obtain a system whose 
classical equilibria are those configurations with the particles lying at relative rest on the vertices of a 
regular tetrahedron of side a. The inertia tensor of the equilibrium configuration is now met 2 diag(l, 1, 1). 
The constants T aa i — e a jiZ a j, a = 1, 2, 3, defining the Eckart gauge are then normalized to !SHf 2 3 = ma 2 . 

This system, unlike that of the previous section, is fully three-dimensional. With e = ^Jfij (mua 2 ), the 
C(e 2 ) perturbation Hi does not commute with the zeroth-order quadratic Hamiltonian, so the 0(e 2 ) correc- 
tions to the unperturbed energies must be computed perturbatively. Diagonalizing TL 1 within eigenspaces 

of Tio is best done numerically, due to the large accidental degeneracies of the unperturbed levels beyond 
the ground state. For that reason, we will restrict ourselves to making only some remarks on the form of 
the 0(e 2 ) Hamiltonian. 

We have six vibrational normal modes, with normal coordinates {SQ a }a=i whose expressions in terms 
of the position vectors {SR a } we omit for brevity. The unperturbed Hamiltonian is given by (JTUIl), with 
/3 now running up to 4, a up to 9, and with a\ — a\ = 1, a 2 = a 2 = a\ = 2 and a 2 = 4. We have 
A^J 1 = 1 / (ma 2 )8ij +C(e 3 ), and the quantum potentials Vi,a and the anharmonic corrections to V( 2 ) starting 
at C(e 3 ) up to constant terms. Thus, to 0(e 2 ) the perturbation Hamiltonian is given by the second term 
in with the momentum operators of 184|) . Dropping constant terms, the C(e 2 ) perturbation can be 

arranged in the form, 

Wx = ^ (s + Af + (e») , A, = \ ± -^Oc^ • (107) 

c— 4 

The normal coordinate we call SQg corresponds to a vibrational mode Tg with Tg a i oc Z a i, i.e., a dilatation 
mode. Explicit computation shows that 5Qgi = 0, and then [A, SQg] — 0. We obtain also [Aj, Y^l=4 $Qa} = 0' 
but [ Aj, V( 2)] 7^ 0. The operator A is not an angular momentum operator, as expected on general grounds 
from (|28bp> . but it turns out to be proportional to one, [2Aj,2Aj] = ieyfe2Afe. We cannot give at present 
necessary and sufficient conditions a many-body system and a rotating frame must satisfy for A, or a multiple 
thereof, to be an angular momentum operator. 



6 Final remarks 

The gauge-invariant approach presented here leads to a general and systematic framework for the quan- 
tization of many-body systems in rotating frames. Our approach naturally incorporates the notions of 
time-dependent symmetry transformations (i.e., gauge transformations), body- frame time-derivatives (co- 
variant derivatives), moving reference frames defined as functions of the particle positions (gauge conditions), 
and of reference-frame singularities (Gribov ambiguities) in a most economical way. It is not, therefore, a 



22 



superfluous formal structure imposed on the physics. The amount of formalism that has been introduced is 
in fact minimal. Rather, we put all those notions within a consistent mathematical framework. 

We have shown that the rotational symmetry of an iV-body system is a gauge symmetry, if we restrict 
ourselves to a fixed angular momentum sector (eq. ©)■ Using gauge invariance we formulated both the 
classical and quantum theory (in the operator approach) in rotating frames defined by linear gauge con- 
ditions. In particular, we explicitly obtained the Hamiltonian operator Ij42(l in terms of position vectors 
referred to the rotating frame, therefore constrained by the gauge conditions. We also showed that the 
orientational degrees of freedom can be eliminated from the formalism, and computed the Hamiltonian op- 
erator and inner product within the corresponding reduced Hilbert space (eqs. <|46[) and (|56[) . resp.). A 
further simplification is obtained by eliminating the Jacobian from the kinetic energy and the inner product, 
leading to the form l|tjtj|) for the Hamiltonian, including the quantum potentials Vi 2 of l|t)5|) . and (|68|l for 
the inner product. The Hamiltonian i|tjtj|) . being Weyl-ordered, can be associated to a generating functional 
in the path-integral approach with mid-point discretization. Those results were extended in section 0] to the 
translation-invariant case, where the system is further reduced by eliminating the center-of-mass degrees of 
freedom. The particular case of quasi- rigid systems was discussed in section [SJ 

The results given in the foregoing apply to a very general class of models comprising all A-particle sys- 
tems in three-dimensional Euclidean space with rotation-invariant potentials. The fact that we computed the 
Hamiltonian operator in terms of position vectors referred to a rotating frame amounts to a purely conven- 
tional choice of coordinates. Once the Hamiltonian has been given in those coordinates it is straightforward 
to transform it to any other coordinate set, as done in sects . and I5"2l there being no need to compute it 
again. For simplicity, however, we restricted ourselves to systems with spin-independent interactions. The 
extension of the formalism to include dynamical spin degrees of freedom should in principle be straightfor- 
ward. We note, in this respect, that in order to obtain half-integer values for the total angular momentum 
of the A-particle system, when appropriate, we should substitute the rigid rotator in by a "rotator" 
with half-integer angular momentum. We did not consider, either, those cases in which it is not possible, or 
desirable, to impose three gauge conditions depending only on particle coordinates. Among those are, e.g., 
one-particle systems (including translation-invariant two-particle systems), for which the gauge conditions 
Ry = Rz = = £,x lead to a description in spherical coordinates, and the conditions Ry = = £,x = £,y 
to cylindrical coordinates. The case of reference frames defined by gauge conditions of a more general form 
than those considered in the previous sections can also be treated by the methods discussed in this paper. 
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A The Laplacian in configuration space 

The kinetic energy of the system of particles considered in section [3] is proportional to the Laplacian in 
configuration space = Yla=i ^ l^'&i w ^ n * ne generalized coordinates {qa\ a =i defined in (|3*T1) . In this 
appendix we compute the expression for in curvilinear coordinates {Q a }a=i defined as follows. For 
1 < a < 3, Q a = 6 a , with 9 a parametrizing U in (fL2")l . and {Q a }a=4 defined by (|3"2"|) . The relation between 
the two sets of coordinates is given implicitly by (|31|l and (|32J) . is then given by the standard expression, 




(A.l) 



dq c dq c 
dQ a dQ b ' 



3iV 



dQ a dQ b 
dq c dq c ' 



(H) = 



M ab = 



J = det 



1/2 



a=l 



In order to obtain an explicit expression for 



we have to build the matrix M, 



r-l 

ab ' 
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For 1 < a < 3 we have dQ a /dq c = d9 a /dq c — J2a=i(d r aj I f ®Qc)(dO a / 'dr a j). From lfT3|) we have, 



ml 



d9 a _ 1 A _! 



U k i = -A i ^^Q. ib 1 m a YbakU] 



(A.2) 



6=1 



where in the second equality we used (|16J) . Thus, 

JV 3 



For 4 < a < 3A we have, from 



a=l b=l 



(A.3) 



N 3 
a - 2^ to2 1 Q /3 J Q r . ~ m2 LamU V Z^ UalU id yy 2 1 dakU k 3 



p 5K 2 Qpi dr aj ft 2 

3=1 J 



(A.4) 



d=l 



where the last equality follows directly from l|15|l and (|16() . Notice that, by definition, O^ 1 is a 3 x 3 matrix 
inverse to Q a i with 1 < a < 3, but in general Q a i£l~l ^ <5 a( j if 4 < a < 37V like in the last term in ljA.41) . 
With ljA.4|) and we obtain, 



in N f N 3 r 

»Va _ \ ~* "^-a ^ rr 1 caj \ " \ " ,-1-1 m a ^ 1 cnj 



a=l d=l 



9^ 



(A.5) 



The matrix elements M~ h x can now be computed, starting with their definition (|A.lfl . and using the orthog- 
onality and completeness relations H3Q(I . and the definition of a ;. 
For 1 < a, b < 3 we get 



3JV 



A C = E = ^ V E ^-n- - a-a-a/- 



(A.6) 



d=l 



with A/" 1 defined by the last equality (compare (|34|) '). For 1 < a < 3, 4 < b < 3 A we get 



(A.7) 



The case 4 < a < 3 A, 1 < b < 3 follows from (|A~7)l by the symmetry of M~ b l . Finally, for 4 < a, 6 < 3A 

3N 



Ko 1 = E 



dQ a dQ b 1 



^ Sab + ^jQaiQbjjVii 1 



dq c dq c m 2 $n 4 



(A. 



M ah x is given by (|A.6|1 - HA.8|) in four blocks 



M- 1 = 



3x3 


3x(3JV-3) 


(3A r -3)x3 


(3Af-3)x(37V-3) 



A~ lt J\f~ 1 A~ 1 


A^Af^tf/m 2 * 


OAf~ 1 A~ 1 /W 2 


l/^ 2 + £3A/'" 1 £3 t /^ 4 



A" 1 * 





o/m 2 





A/"" 1 








Af- 1 /^ 2 



A- 1 










(A.9) 



From the last equality we obtain J = m 3 \A\J, with |A| = det(A) and J = (det(A/")) 1/2 . Substituting fO 
and J into (|A. If) we obtain V 2 = -2(Hn - V), with Hjv given by l|33jl. 
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A.l Reduced Laplacian and Weyl ordering 

In order to eliminate the factors of J from V 2 q in (| A. 1|) . and to Weyl-order it, we write, 

^ -= J1, ^ J - 1 ' 2 - £ fa'a&s + " ^ £ (4: ( M -0) • < A - 10) 

a, o—l a, o—l 

Defining the Weyl-ordered differential operator, 

' ah dQ a dQ b ) w ~ A ab dQ a dQ b + AdQ a ab dQ b + 4 8Q b ab dQ a + 4 dQ a dQ b ab { ' ' 



(|A.10p becomes 



9 I ao dQadQb ) 4 4^ 5Q a 9Q ^ WQa V ^ 9Q h 



a,b=l x " a,6=l ~ ~ v "~ a,b=l 



(A.12) 



The last two terms on the r.h.s. are multiplicative operators which, up to a constant factor, constitute the 
quantum potential. The last one can be considerably simplified by using the second line of (|A.1(I to write, 



1 ( 9 ( M - l9J1/2 \\ - 1 V J_J_(^Lji/2 

oil Jl/2 W« V a " dQ b J J 2 a f^ =i Ji/2 dQ a \ dq c \dQ b dq c 



1 / d 8Q a \ ( d dQ b \ , 19Q a / d 2 dQ 



3JV 

E 

',,fe,C— 1 



I dQ a ( d dQ b \ ^ 8Q d ( d dq, 



2 \dQ a dq c J \dQ b dq c J 2 dq c \dQ a dQ b dq c 

E 



(A.13) 



4 <9g c \dQ b dq c J f-r 1 , dq c > \dQ a dQ d 

c ,e£— 1 

where for the last equality we used the analog of (|35|1 . 



c,a— 1 

The last term in (|A.13|I can be simplified using 

E dQd ( d dq c > \ _ y f d dQ d \ dq 
Fin., \ f)0. Fin, I ~ Z^ 



. dq c > \dQadQ d J f-f \dQa dq c > ) dQd 

d—1 c'd—1 



V 9Q a ( d dQ d \ ^dQ a ( d dQ d 



and 



to obtain, 



v 4 \dQ a dq c J \dQ b dq c J 2 dq c \dQ a dQ b dq, 

a,o,c— 1 

The second term on the r.h.s. of (|A.12Jl can be rewritten as, 

1 ^ 9 2 M- b x 1 ^ / / d 2 dQ a \ dQ b f d dQ a \ ( d dQ 



Z-j dq c V dQa dq d ) Zv dq d \ dQ a dq 

_ ^ (I ( d dQ a \ ( d dQ b \ IdQa ( d 2 dQ b ^ v 
~ Z^ 7 [on ST ^ + 9 ^ ^n, ^ ' lA - i0j 



4 an an, a Z^ 



4 4^ dQ a dQ b 4 4^ . V \9Q a dQ b dq c J dq c \dQ b dq c J \dQ a 

a, o—l a,o,c— 1 
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/ d dQ a \ ( d dQ 



\dQ a dq c ) \dQ b dq c 
Thus, finally, substituting (|A. 15|) and JA.16|) into <|A.12() . 

3N , a a \ , 3W 

V «- 2. \ M a b g QadQt 



,.fc=l 



11' 



3N 



,.fc.c=l 



9 5Q a \ / 9 dQt 



dQb dq c ) \8Q a 



A. 2 The Quantum potential 

The kinetic energy in l|58|) is — 1/2V 2 ,. Thus, from (|A. 1T|> we have, 



:SiY 



Vq 



- E 

8 f-^ 

a,b,c— 1 



/ 9 dQ a \ f d dQ 



\dQ b dq c J \dQa dq 



(A.16) 



(A.17) 



(A.18) 



The evaluation of Vq is straightforward, though rather laborious. We closely follow the analogous computa- 
tion of UJ, which leads to a compact expression for Vq. In order to compute Vq we split it into three terms, 



^ 3 3N 3N , q 8Q b \ ( d 86 \ 1 3 3Ar 

= 4^S? v^^fJ fe^rJ ' Va> = 8 5§ 



(A.19) 



a=l b=4 c=l x " ■"■/ \ J<-/ a,b=l c=l 

which we consider separately. 
A.2.1 Vq 

From l|A.5|> we get dQb/dq c and thence, 

IV „ 3 



89 a dq c J V d9 b 



f) - - g £ g (£0^) . (a. 2 o, 

Substituting this expression, and the corresponding one with a, 6 interchanged, into Vq and using the 
completeness and orthogonality relations (|30|l . and (|40|l . we get, 

37V n 3 2 

V Q° =lE EE r< ^ r ^' ) (P 7 ^0 om 0- x d ) . (A.21) 

es,6=4 (8,7=1 <i=l 

In (|A.21|) we can extend the summations over a and 6 down to 1, due to 1)24(1. Expanding the definitions @ 
of £lbg and £J am and using completeness, (|3L)|) . 

^ iV 3 

^Qo = "g E E^^' ffs ^m P (PptRyQ^) (P^RppQ^) . (A.22) 

/3,7=ld=l 

It is not difficult to check that if in this equation we expand the expression (I29f) for Ppi and P 7 i', the 
contribution due to the second term in (|29|1 vanishes, and we get, 

N 3 

^ = sEE (afe^ 1 ) (4^°^) ■ (A - 23) 



2G 



Using the definition J5J of Q. a i, the derivatives can be evaluated to give, 

id^i R ' fsQ '^) = ( S ^ SslS <> k ~ R ~< s E ^gam^ a0n e nk A H^ 1 (A.24) 
and similarly (8/dRyi'RppQ.^). Thus, 

V Qo = g E ( X! ^d^fcd^M ) £ l'gs £ lmp ( S 01 5 sl 8 gk - R JS Qga™,f3T a p n £ nk t j 
(8,7=1 \<i=l / V o=l / 

. (A.25) 



b=l 



Identifying the first parenthesis in this expression with M k} ^ as defined in i|34|) . it is straightforward to 
rearrange the factors to obtain Vq = V2, with V2 given in (|65|) . 

A.2.2 V Ql and Vq 2 

With the derivatives d9 a /dq c given in l|A.3|) . using the completeness relation l|30[) we easily get, 

3 iv 3 / o \ / q \ 

^ = gEE E "»«r«i 1 a*r lfa a P 0^£2^ ^-A- 6 V«) ^Va 1 ^ • (A.26) 

a,6=la=ldi,d 2 =l ^ a / \ & / 

We can write the derivatives dll/d9 a in terms of A a i using l|13|) to obtain the expanded form, 



v 1 V (rtft)-i{}A aA * aA £ 

Q2 ~*h Ih ld qd m « & 

a,o— 1 \a— 1 / 



1 3 N 3 ( dk~ l (9A _1 \ 

s E E E ^r.^r^n^o-^ -g^K 1 ^ - -Q^Kb A ™ (A.27) 



a, 6—1 a— 1 rfi,d 2 — 1 
iV 3 



o E/ E/ m ard 1 Q/ c rd 2 api3; [ i i n 9( i 2 £p;r£fegr 



8 

Q=l di ,d2 = l 



The second line can be evaluated from the commutators (|19fl . [Li,Lj] = —iSij k L k . Substituting into these 
commutators the expression (1221) for Li, we are led to 



d=l 

Thus, 1A.27II can be rewritten as, 



1 3 dA~ 1 QA 1 1 N 3 

^2 = g E ^ 1 ~de^~~d0Y + 8 ^ ^ maTd ^ anS:i qd2 ( £ nlr£kqr + £qll'£kl'n) Qf/T^ak ■ (A.29) 
a : b— 1 a— ldi,d 2 — 1 

The factor in parenthesis in the second term equals £ n qi'£kU', but we will refrain from simplifying it. After 
renaming dummy indices we can rewrite l|A.29l) as, 

V Q? = n E ^U' 1 J) 11 ' -p^T- + o E m « E {&l>l 2 T d 2aq £qln - £UlqQqd 2 T d2an) Q^/j 'd ia k£ kV n , 



d9 a d9 a , 

a.a' — l 



(A.30) 



27 



which is exactly analogous to (6.8) of |SJ. 

With the derivatives (|A. 3|) and (|A.5|I . and using completeness ()3(J[1 . we can write 

1 3 3N n / 3 \ <9[7 3 SCT 1 

= _ i H EI 51 TrI ( r6ar ~ EI &b q Qglr d2ar 1 ^ T dian A-^U ni -^- . (A.31) 

a=l f)=4 Q=l \ d 2 = l / a di = l 

The derivative dU r i/d9 a can be written in terms of A Q i with (|13() . On the other hand, using (|39|l to write 
dQ.Ji^ /dQb in terms of dQ^^/dRpi and evaluating the latter from we arrive at 

1 N ( 3 \ 3N m m f 3 \ 

^ = -- ^ EI n^r^pSp,™ (n^^Td^n) ^ 2 I r fcQ r - X! £6«£gd 2 r <w 1 £ri ,„ . 

a, 0=1 \c,d 1 = l J b=4 V d 2 = l ) 

(A.32) 

Expanding 0b g and using completeness to evaluate the sum over b, with T>f q defined in (|59|l we get 
I N 3 , / 3 \ 

VQi = - J X! m « X! (^J^crc/Jp^p/m) (^mdi^icm) I ^a/^ir ~ Ej 771 P V Iq&qd^ d 2 ar j £rf'n , (A.33) 

an expression which is exactly analogous to eq. (6.10) of [jy. Combining this last expression for Vqi and 
(|A.29|I for Vq2, we can write, 

1 3 dA~ 1 <9A 1 1 ^ 3 

VQl + V Q2 = g X! -^r 1 dg a ' ~Q0 2 7 ~ g X! mQ X! 5 3i 7 c r ca;^ m d r 'iQ™( £ ™i'p e p/m + £nZp£p;'m) 

a, a' — 1 a a a— 1 c.d— 1 . . , 

(A.34 

^ AT 3 3 v ; 

"7 E m a m /3 EI iCl i"c r ^P £ pim £3 mo!i rd i a '" £: '"''n EI V lq^qd\ V d 2 an ■ 



4 

a,/3=l c,di = l d 2 = l 



The summand on the second line of this equation is the product of (e p / m I?^) times an expression antisym- 
metric in m and q. Thus, we can replace (e p / m I?^) — > l/2(£ p j m D^ — s p i q T>[ m ) — l/2s q i m 'Df p , and the second 
line of l)A.34|l becomes, 

I f N 3 \ N 3 

~g E^E^'c 11 *^ £ ?I™E m ° EI ^md^di^rSW'n^^r^an ■ (A.35) 
\3=1 c=l y a=l rfi,d 2 =l 

From the definitions (j59(l and (0, we see that the factor within parentheses in l|A.35|l reduces to (—5w). 
Making that simplification in (|A.35|) and substituting the result for the second line of (|A.34|I we finally get, 

a,a' — 1 



with Vi given by JSHJl. 



B AfRne transformations and covariant derivatives 

The transformations (|69|l do not depend on r a linearly but affinely, i.e., {c\r a + C2rp)' — cir a ' + C2rp iff 
Ci + c% = 1. Covariant derivatives are defined so that they transform under time-dependent transformations 
in the same way as ordinary derivatives transform under time-indepent transformations. Thus, we define 
D t r a as in (fTTTf) . so that (D t r a )' = UD t r a , but D t (c\r a + C2Tp) = ciD t r a + C2D t r p iff ci + c 2 = 1. Similarly, 
the rule for the derivative of a vector product is not the usual one, D t (r a Arp) — (D t r a ) Ar^ + r Q A (D t rp) + 
(fa ~ r p) A p — p. Since D t r a transforms linearly under gauge transformations, we define, 

D t D t r a = ^-(D t r a ) - £(D t r a ) , (D t D t r a )' = UD t D t r a . (B.l) 
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With l|7U|) and (jB.lfl . the eqs. of motion from Lagrangian (|72[1 take the form m a D t D t r a = — V Q V as in QJ. 

With the definition (jSJ) for the angular momentum I and the transformation law Q69|l we have I' = 
Ul + Mu A (UD t r cm ), with r cm the center-of-mass position vector. Thus, we define, 

D t l = i-£l-MpA{D t r cal ) , (B.2) 

so that (D t l)' = UD t l + MuA (UD t D t r cm ) with D t D t r cm defined as in [JB.lfl and {7UJ. The center-of-mass 
angular momentum i cm = Afr cm A D t r cra transforms in the same way as Z, and its covariant derivative is 
defined as in l|B.2(l . From the eqs. of motion for r a we then get, 

D t l = , D t D t r cm = , D t l cm = . (B.3) 

Therefore, D t l — D t l C m — (d/dt — £)(Z — Z cm ) = 0, which is (|73|) and which, together with the antisymmetry 
of immediately leads to d/dt(l — i cm ) 2 = 0. Furthermore, (I — i C m)' = U(l — Z C m) so that (I — l cm ) 2 is 
invariant under gauge transformations, i.e., frame-independent. 
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